<!DOCTYPE html><html lang="zh-CN" data-theme="light"><head><meta charset="UTF-8"><meta http-equiv="X-UA-Compatible" content="IE=edge"><meta name="viewport" content="width=device-width, initial-scale=1.0, maximum-scale=1.0, user-scalable=no"><title>基于Fortran90的Rpstar3大气辐射传输系统源代码文档 | 洛沐の人间客栈</title><meta name="keywords" content="辐射传输"><meta name="author" content="洛沐,guojiaxiang0820@gmail.com"><meta name="copyright" content="洛沐"><meta name="format-detection" content="telephone=no"><meta name="theme-color" content="#ffffff"><meta name="description" content="这是日本东京大学辐射传输源码官方文档的中文版翻译，仅做学习之用。">
<meta property="og:type" content="article">
<meta property="og:title" content="基于Fortran90的Rpstar3大气辐射传输系统源代码文档">
<meta property="og:url" content="https://www.guojxblog.cn/archives/8c1f6a98.html">
<meta property="og:site_name" content="洛沐の人间客栈">
<meta property="og:description" content="这是日本东京大学辐射传输源码官方文档的中文版翻译，仅做学习之用。">
<meta property="og:locale" content="zh_CN">
<meta property="og:image" content="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/%E5%9F%BA%E4%BA%8EFortran90%E7%9A%84Rpstar3%E8%BE%90%E5%B0%84%E4%BC%A0%E8%BE%93%E6%A8%A1%E5%9E%8B%E6%BA%90%E4%BB%A3%E7%A0%81%E6%96%87%E6%A1%A3/rpstar.jpeg">
<meta property="article:published_time" content="2022-05-21T03:05:40.000Z">
<meta property="article:modified_time" content="2023-01-18T14:14:11.756Z">
<meta property="article:author" content="洛沐">
<meta property="article:tag" content="辐射传输">
<meta name="twitter:card" content="summary">
<meta name="twitter:image" content="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/%E5%9F%BA%E4%BA%8EFortran90%E7%9A%84Rpstar3%E8%BE%90%E5%B0%84%E4%BC%A0%E8%BE%93%E6%A8%A1%E5%9E%8B%E6%BA%90%E4%BB%A3%E7%A0%81%E6%96%87%E6%A1%A3/rpstar.jpeg"><link rel="shortcut icon" href="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/RS_icon.jpg"><link rel="canonical" href="https://www.guojxblog.cn/archives/8c1f6a98"><link rel="preconnect" href="//cdn.jsdelivr.net"/><link rel="preconnect" href="//busuanzi.ibruce.info"/><meta name="baidu-site-verification" content="code-zGPrfnFp7y"/><link rel="stylesheet" href="/css/index.css"><link rel="stylesheet" href="https://cdn.jsdelivr.net/npm/@fortawesome/fontawesome-free@6/css/all.min.css" media="print" onload="this.media='all'"><link rel="stylesheet" href="https://cdn.jsdelivr.net/npm/node-snackbar/dist/snackbar.min.css" media="print" onload="this.media='all'"><link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/animate.css/4.1.1/animate.min.css" media="print" onload="this.media='all'"><script>const GLOBAL_CONFIG = { 
  root: '/',
  algolia: undefined,
  localSearch: {"path":"search.xml","languages":{"hits_empty":"找不到您查询的内容：${query}"}},
  translate: {"defaultEncoding":2,"translateDelay":0,"msgToTraditionalChinese":"繁","msgToSimplifiedChinese":"簡"},
  noticeOutdate: {"limitDay":180,"position":"top","messagePrev":"此文中创建于","messageNext":"天之前，请以最新文章为准！"},
  highlight: {"plugin":"highlighjs","highlightCopy":true,"highlightLang":true,"highlightHeightLimit":500},
  copy: {
    success: '复制成功',
    error: '复制错误',
    noSupport: '浏览器不支持'
  },
  relativeDate: {
    homepage: false,
    post: true
  },
  runtime: '',
  date_suffix: {
    just: '刚刚',
    min: '分钟前',
    hour: '小时前',
    day: '天前',
    month: '个月前'
  },
  copyright: {"limitCount":50,"languages":{"author":"作者: 洛沐","link":"链接: ","source":"来源: 洛沐の人间客栈","info":"著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。"}},
  lightbox: 'mediumZoom',
  Snackbar: {"chs_to_cht":"你已切换为繁体","cht_to_chs":"你已切换为简体","day_to_night":"你已切换为深色模式","night_to_day":"你已切换为浅色模式","bgLight":"#49b1f5","bgDark":"#1f1f1f","position":"bottom-left"},
  source: {
    justifiedGallery: {
      js: 'https://cdn.jsdelivr.net/npm/flickr-justified-gallery@2/dist/fjGallery.min.js',
      css: 'https://cdn.jsdelivr.net/npm/flickr-justified-gallery@2/dist/fjGallery.min.css'
    }
  },
  isPhotoFigcaption: true,
  islazyload: false,
  isAnchor: true
}</script><script id="config-diff">var GLOBAL_CONFIG_SITE = {
  title: '基于Fortran90的Rpstar3大气辐射传输系统源代码文档',
  isPost: true,
  isHome: false,
  isHighlightShrink: false,
  isToc: true,
  postUpdate: '2023-01-18 22:14:11'
}</script><noscript><style type="text/css">
  #nav {
    opacity: 1
  }
  .justified-gallery img {
    opacity: 1
  }

  #recent-posts time,
  #post-meta time {
    display: inline !important
  }
</style></noscript><script>(win=>{
    win.saveToLocal = {
      set: function setWithExpiry(key, value, ttl) {
        if (ttl === 0) return
        const now = new Date()
        const expiryDay = ttl * 86400000
        const item = {
          value: value,
          expiry: now.getTime() + expiryDay,
        }
        localStorage.setItem(key, JSON.stringify(item))
      },

      get: function getWithExpiry(key) {
        const itemStr = localStorage.getItem(key)

        if (!itemStr) {
          return undefined
        }
        const item = JSON.parse(itemStr)
        const now = new Date()

        if (now.getTime() > item.expiry) {
          localStorage.removeItem(key)
          return undefined
        }
        return item.value
      }
    }
  
    win.getScript = url => new Promise((resolve, reject) => {
      const script = document.createElement('script')
      script.src = url
      script.async = true
      script.onerror = reject
      script.onload = script.onreadystatechange = function() {
        const loadState = this.readyState
        if (loadState && loadState !== 'loaded' && loadState !== 'complete') return
        script.onload = script.onreadystatechange = null
        resolve()
      }
      document.head.appendChild(script)
    })
  
      win.activateDarkMode = function () {
        document.documentElement.setAttribute('data-theme', 'dark')
        if (document.querySelector('meta[name="theme-color"]') !== null) {
          document.querySelector('meta[name="theme-color"]').setAttribute('content', '#0d0d0d')
        }
      }
      win.activateLightMode = function () {
        document.documentElement.setAttribute('data-theme', 'light')
        if (document.querySelector('meta[name="theme-color"]') !== null) {
          document.querySelector('meta[name="theme-color"]').setAttribute('content', '#ffffff')
        }
      }
      const t = saveToLocal.get('theme')
    
          const now = new Date()
          const hour = now.getHours()
          const isNight = hour <= 6 || hour >= 18
          if (t === undefined) isNight ? activateDarkMode() : activateLightMode()
          else if (t === 'light') activateLightMode()
          else activateDarkMode()
        
      const asideStatus = saveToLocal.get('aside-status')
      if (asideStatus !== undefined) {
        if (asideStatus === 'hide') {
          document.documentElement.classList.add('hide-aside')
        } else {
          document.documentElement.classList.remove('hide-aside')
        }
      }
    
    const detectApple = () => {
      if(/iPad|iPhone|iPod|Macintosh/.test(navigator.userAgent)){
        document.documentElement.classList.add('apple')
      }
    }
    detectApple()
    })(window)</script><link rel="stylesheet" href="/css/style.css"><link rel="stylesheet" href="/css/iconfont.css"><link rel="stylesheet" href="/css/font-awesome.min.css"><link rel="stylesheet" href="https://cdn.jsdelivr.net/gh/sviptzk/StaticFile_HEXO@latest/butterfly/css/macblack.css"><script src="/live2d-widget/autoload.js"></script><meta name="generator" content="Hexo 6.1.0"></head><body><div id="web_bg"></div><div id="sidebar"><div id="menu-mask"></div><div id="sidebar-menus"><div class="avatar-img is-center"><img src="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/loading3.gif" data-original="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/touxiang.jpeg" onerror="onerror=null;src='/img/friend_404.gif'" alt="avatar"/></div><div class="site-data is-center"><div class="data-item"><a href="/archives/"><div class="headline">文章</div><div class="length-num">29</div></a></div><div class="data-item"><a href="/tags/"><div class="headline">标签</div><div class="length-num">18</div></a></div><div class="data-item"><a href="/categories/"><div class="headline">分类</div><div class="length-num">7</div></a></div></div><hr/><div class="menus_items"><div class="menus_item"><a class="site-page" href="/"><i class="fa-fw fas fa-home"></i><span> 首页</span></a></div><div class="menus_item"><a class="site-page group hide" href="javascript:void(0);" rel="external nofollow noreferrer"><i class="fa-fw fa fa-graduation-cap"></i><span> 文章</span><i class="fas fa-chevron-down"></i></a><ul class="menus_item_child"><li><a class="site-page child" href="/archives/"><i class="fa-fw fa fa-folder-open"></i><span> 归档</span></a></li><li><a class="site-page child" href="/categories/"><i class="fa-fw fa fa-archive"></i><span> 分类</span></a></li><li><a class="site-page child" href="/tags/"><i class="fa-fw fa fa-tags"></i><span> 标签</span></a></li></ul></div><div class="menus_item"><a class="site-page group hide" href="javascript:void(0);" rel="external nofollow noreferrer"><i class="fa-fw fa-solid fa-heart-pulse"></i><span> 生活</span><i class="fas fa-chevron-down"></i></a><ul class="menus_item_child"><li><a class="site-page child" href="/books/"><i class="fa-fw fas fa-book-reader"></i><span> 读书</span></a></li><li><a class="site-page child" href="/photos/"><i class="fa-fw fa fa-camera-retro"></i><span> 相册</span></a></li><li><a class="site-page child" href="/music/"><i class="fa-fw fa fa-music"></i><span> 音乐</span></a></li><li><a class="site-page child" href="/movies/"><i class="fa-fw fas fa-video"></i><span> 影视</span></a></li></ul></div><div class="menus_item"><a class="site-page group hide" href="javascript:void(0);" rel="external nofollow noreferrer"><i class="fa-fw fa-solid fa-keyboard"></i><span> 动态</span><i class="fas fa-chevron-down"></i></a><ul class="menus_item_child"><li><a class="site-page child" href="/comment/"><i class="fa-fw fa fa-paper-plane"></i><span> 留言</span></a></li><li><a class="site-page child" href="/notes/"><i class="fa-fw fas fa-feather-alt"></i><span> 随笔</span></a></li><li><a class="site-page child" href="/ideas/"><i class="fa-fw fa-solid fa-lightbulb"></i><span> 妙想</span></a></li></ul></div><div class="menus_item"><a class="site-page" href="/link/"><i class="fa-fw fa fa-link"></i><span> 友链</span></a></div><div class="menus_item"><a class="site-page" href="/about/"><i class="fa-fw fas fa-heart"></i><span> 关于</span></a></div></div></div></div><div class="post" id="body-wrap"><header class="post-bg" id="page-header" style="background-image: url('https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/%E5%9F%BA%E4%BA%8EFortran90%E7%9A%84Rpstar3%E8%BE%90%E5%B0%84%E4%BC%A0%E8%BE%93%E6%A8%A1%E5%9E%8B%E6%BA%90%E4%BB%A3%E7%A0%81%E6%96%87%E6%A1%A3/rpstar.jpeg')"><nav id="nav"><span id="blog_name"><a id="site-name" href="/">洛沐の人间客栈</a></span><div id="menus"><div id="search-button"><a class="site-page social-icon search"><i class="fas fa-search fa-fw"></i><span> 搜索</span></a></div><div class="menus_items"><div class="menus_item"><a class="site-page" href="/"><i class="fa-fw fas fa-home"></i><span> 首页</span></a></div><div class="menus_item"><a class="site-page group hide" href="javascript:void(0);" rel="external nofollow noreferrer"><i class="fa-fw fa fa-graduation-cap"></i><span> 文章</span><i class="fas fa-chevron-down"></i></a><ul class="menus_item_child"><li><a class="site-page child" href="/archives/"><i class="fa-fw fa fa-folder-open"></i><span> 归档</span></a></li><li><a class="site-page child" href="/categories/"><i class="fa-fw fa fa-archive"></i><span> 分类</span></a></li><li><a class="site-page child" href="/tags/"><i class="fa-fw fa fa-tags"></i><span> 标签</span></a></li></ul></div><div class="menus_item"><a class="site-page group hide" href="javascript:void(0);" rel="external nofollow noreferrer"><i class="fa-fw fa-solid fa-heart-pulse"></i><span> 生活</span><i class="fas fa-chevron-down"></i></a><ul class="menus_item_child"><li><a class="site-page child" href="/books/"><i class="fa-fw fas fa-book-reader"></i><span> 读书</span></a></li><li><a class="site-page child" href="/photos/"><i class="fa-fw fa fa-camera-retro"></i><span> 相册</span></a></li><li><a class="site-page child" href="/music/"><i class="fa-fw fa fa-music"></i><span> 音乐</span></a></li><li><a class="site-page child" href="/movies/"><i class="fa-fw fas fa-video"></i><span> 影视</span></a></li></ul></div><div class="menus_item"><a class="site-page group hide" href="javascript:void(0);" rel="external nofollow noreferrer"><i class="fa-fw fa-solid fa-keyboard"></i><span> 动态</span><i class="fas fa-chevron-down"></i></a><ul class="menus_item_child"><li><a class="site-page child" href="/comment/"><i class="fa-fw fa fa-paper-plane"></i><span> 留言</span></a></li><li><a class="site-page child" href="/notes/"><i class="fa-fw fas fa-feather-alt"></i><span> 随笔</span></a></li><li><a class="site-page child" href="/ideas/"><i class="fa-fw fa-solid fa-lightbulb"></i><span> 妙想</span></a></li></ul></div><div class="menus_item"><a class="site-page" href="/link/"><i class="fa-fw fa fa-link"></i><span> 友链</span></a></div><div class="menus_item"><a class="site-page" href="/about/"><i class="fa-fw fas fa-heart"></i><span> 关于</span></a></div></div><div id="toggle-menu"><a class="site-page"><i class="fas fa-bars fa-fw"></i></a></div></div></nav><div id="post-info"><h1 class="post-title">基于Fortran90的Rpstar3大气辐射传输系统源代码文档</h1><div id="post-meta"><div class="meta-firstline"><span class="post-meta-date"><i class="far fa-calendar-alt fa-fw post-meta-icon"></i><span class="post-meta-label">发表于</span><time class="post-meta-date-created" datetime="2022-05-21T03:05:40.000Z" title="发表于 2022-05-21 11:05:40">2022-05-21</time><span class="post-meta-separator">|</span><i class="fas fa-history fa-fw post-meta-icon"></i><span class="post-meta-label">更新于</span><time class="post-meta-date-updated" datetime="2023-01-18T14:14:11.756Z" title="更新于 2023-01-18 22:14:11">2023-01-18</time></span><span class="post-meta-categories"><span class="post-meta-separator">|</span><i class="fas fa-inbox fa-fw post-meta-icon"></i><a class="post-meta-categories" href="/categories/%E7%A8%8B%E5%BA%8F%E4%BB%A3%E7%A0%81/">程序代码</a></span></div><div class="meta-secondline"><span class="post-meta-separator">|</span><span class="post-meta-wordcount"><i class="far fa-file-word fa-fw post-meta-icon"></i><span class="post-meta-label">字数总计:</span><span class="word-count">7.2k</span><span class="post-meta-separator">|</span><i class="far fa-clock fa-fw post-meta-icon"></i><span class="post-meta-label">阅读时长:</span><span>28分钟</span></span><span class="post-meta-separator">|</span><span class="post-meta-pv-cv" id="" data-flag-title="基于Fortran90的Rpstar3大气辐射传输系统源代码文档"><i class="far fa-eye fa-fw post-meta-icon"></i><span class="post-meta-label">阅读量:</span><span id="busuanzi_value_page_pv"></span></span></div></div></div></header><main class="layout" id="content-inner"><div id="post"><article class="post-content" id="article-container"><blockquote>
<p>本文程序在Rstar7与Pstar3的基础上于2020. 3.31 创建，并更新于2020. 8. 7。 英文文件：<a target="_blank" rel="noopener external nofollow noreferrer" href="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/%E5%9F%BA%E4%BA%8EFortran90%E7%9A%84Rpstar3%E8%BE%90%E5%B0%84%E4%BC%A0%E8%BE%93%E6%A8%A1%E5%9E%8B%E6%BA%90%E4%BB%A3%E7%A0%81%E6%96%87%E6%A1%A3/R_Pstar3_Readme.docx">R-Pstar_Readme</a></p>
<p>主要开发者：Miho SEKIGUCHI、Yoshifumi Ota和Teruyuki NAKAJIMA<br>贡献者列表如下:</p>
</blockquote>
<p><img src="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/loading3.gif" data-original="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/%E5%9F%BA%E4%BA%8EFortran90%E7%9A%84Rpstar3%E8%BE%90%E5%B0%84%E4%BC%A0%E8%BE%93%E6%A8%A1%E5%9E%8B%E6%BA%90%E4%BB%A3%E7%A0%81%E6%96%87%E6%A1%A3/Contributor.png" alt="R-Pstar系统贡献者"></p>
<h1 id="“R-Pstar3”简介？"><a href="#“R-Pstar3”简介？" class="headerlink" title="“R-Pstar3”简介？"></a>“R-Pstar3”简介？</h1><p>“Pstar3”是用于在大气-陆地&#x2F;海洋耦合系统中模拟极化辐射场的通用软件包。”R-pstar3’是一个类似Rstar的并用于合并和集成Rstar和Pstar的Pstar包。该模拟基于离散纵坐标&#x2F;矩阵算子法在平面平行问题中的应用分为若干均匀层的大气，下伏地面或海洋表面。Pstar3还能够处理辐射传输与平面平行海洋耦合（Ota et al.，2010）。</p>
<p>最初的Pstar3最初是为观测温室气体而开发的日本国家环境研究所卫星（GOSAT）项目（NIES），由日本JAXA&#x2F;EORC运营。R-Pstar3包通过<a target="_blank" rel="noopener external nofollow noreferrer" href="http://157.82.240.167/~clastr/">Open CLASTR</a>项目分发。</p>
<h1 id="软件包程序清单"><a href="#软件包程序清单" class="headerlink" title="软件包程序清单"></a>软件包程序清单</h1><blockquote>
<ul>
<li>README: 自述文档。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>Makefile: 生成可执行文件的配置文件，可以利用此文档一键编译。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>Mkinclude: Makefile中配置程序中使用的的编译器设置。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>main.f90: R-pstar3软件包主程序。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li><p>data: R-pstar3输入参数样本文件（默认值：指向data0的符号链接）。</p>
<blockquote>
<ul>
<li>data0: 朗伯表面的非偏振辐射</li>
<li>data10: 海洋表面的非偏振辐射（海洋层&#x3D;0）</li>
<li>data1: 海洋表面的非偏振辐射（海洋层&#x3D;2）</li>
<li>data0p:  朗伯表面的偏振辐射</li>
<li>data10p: 海洋表面偏振辐射（海洋层&#x3D;0）</li>
<li>data1p: 海洋表面偏振辐射（海洋层&#x3D;2）</li>
</ul>
</blockquote>
</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li><p>out_test: R-pstar3输出文件的示例，文件命名对应于data文件包：data0,10,1,0p,10p,1p–&gt;out0,10,1,0p,10p,1p。</p>
<p>如果你想得到同样的结果，你应该选择Frohlich和Shaw 1980并在路径LBR&#x2F;pstr4.f90文件中注释掉第528行-第529行（我看到源码中已经注释了）。</p>
</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li><p>这些数据库文件是从ps3r_dbs_200329.tar.gz中提取的。这是一个数据哭的压缩包，需要解压之后将提取的文件放在项目根目录之中。</p>
<blockquote>
<ul>
<li><p>MLATMD: 大气剖面（廓线）文件</p>
</li>
<li><p>AERDB: Rstar7的粒子剖面（廓线）文件</p>
</li>
<li><p>PKRNL.OUT: Mie散射核心文件（增加Mie散射的角度网格）</p>
</li>
<li><p>pyphsf_vis.dat: 六边形固体冰柱数据（近红外和可见光）</p>
</li>
<li><p>pyphsf_ir.dat: 六角形固体冰柱数据（红外）</p>
<p>以上2个文件的字节顺序为little-endian，这些文件与极化辐射不对应。</p>
</li>
<li><p>dkrnl&#x2F;PKRNL.OUT_aspXXX: 球体散射数据</p>
<p>XXX表示宽高比；宽高比&#x3D;X.XX，也就是说如果文件命名为PKRNL.OUT_asp033，那么宽高比就是0.33，这些文件与非极化辐射不对应</p>
</li>
</ul>
</blockquote>
</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>ckd.g.ch_2_1e3_big: big endian机制的k分布文件</li>
<li>ckd.g.ch_2_1e3_ltl: little endian机制的k分发文件</li>
<li>ckd.g.ch_2_2e3_big: big endian机制的精细k分布文件</li>
<li>ckd.g.ch_2_2e3_ltl: little endian机制的精细k分发文件</li>
</ul>
</blockquote>
<h1 id="源码运行与编译"><a href="#源码运行与编译" class="headerlink" title="源码运行与编译"></a>源码运行与编译</h1><p>a、 请仔细阅读本文件。<br>b、 复制计算机目录中的所有文件并将数据库文件放在同一目录中。<br>c、 选择普通或更细的波数栅格，并复制或链接一个文件到新文件’ckd.g.ch_2’。还应指示正常（1e3）或主程序中DINTVL的更精细（2e3）网格。<br>d、 如果要计算球体粒子，选择一个球体散射数据文件，以备最佳使用并在main7.f90中修改文件名。<br>e、 在Mkinclude中修改编译器设置。用户需要指定Fortran编译器（FC）及其选项（FFLAGS和FMATH）。使用FFLAGS进行编译器优化，和用于链接LAPACK库的FMATH（见下文）。</p>
<blockquote>
<p><em>关于FMATH</em><br>Pstar4需要LAPACK和相关例行程序（<em>1）。要构建Pstar4的可执行文件，编译器选项“FMATH”必须为在Mkinclude中指定，以链接LAPACK库。示例选项<br>对于典型的编译器，可以在Mkinclude中找到。强烈建议用户准备优化的LAPACK库用于他们自己的计算机架构（</em>2），并使用“FMATH”链接它选项，因为Pstar4的计算效率在很大程度上取决于LAPACK库，尤其是在IPOL&#x3D;4的情况下。<br>（<em>1）<a target="_blank" rel="noopener external nofollow noreferrer" href="http://www.netlib.org/lapack/">LAPACK — Linear Algebra PACKage</a><br>（</em>2）例如英特尔数学内核库或AMD核心数学库。</p>
<p>注：</p>
</blockquote>
<p>f、 在根目录下输入命令</p>
<figure class="highlight bash"><table><tr><td class="gutter"><pre><span class="line">1</span><br></pre></td><td class="code"><pre><span class="line">make</span><br></pre></td></tr></table></figure>

<p>开始编译，生成Makefile中所要求的执行文件，默认为a.out(MacOS or Linux)或者a.exe(Win)。然后找到此执行文件（在Mkinclude文件中可以设置或者修改执行文件的名称，默认为TARGET &#x3D; a.out）。<br>g、 运行执行文件，输入命令：.&#x2F;的意思是执行当前目录下的执行文件+文件名a.out</p>
 <figure class="highlight bash"><table><tr><td class="gutter"><pre><span class="line">1</span><br></pre></td><td class="code"><pre><span class="line">./a.out</span><br></pre></td></tr></table></figure>

<p>回车之后就开始运行执行文件了。会得到一个输出文件out，其中包括计算结果。输出结果out的内容与out0文件的内容相同，如果有差异可能就是四舍五入误差的容许误差，或者计算机的精度差异。代码正常运行，具有自动双精度选项的双精度运行，fortran编译器可以更好地保持准确性，有时会出现输入&#x2F;输出错误，检查在主程序的最开始处打开语句。<br>h、 对于使用R-Pstar3的应用程序，可以更改数据或者输入参数。改变<strong>大气参数和&#x2F;或粒子光学特性</strong>，必须更改<strong>MLATMD和&#x2F;或AERDB</strong>。请参阅后面的一节。<br>当你计算多案例时，首个案例会非常慢，因为代码需要读取大文件PKRNL，所以制造连续计算，无需拆分为多个作业。<br>i、 您可能需要更改应用程序的参数设置，在这种情况下，应该更改主程序中的模块“paras”。</p>
<h1 id="实验参数文件：数据"><a href="#实验参数文件：数据" class="headerlink" title="实验参数文件：数据"></a>实验参数文件：数据</h1><p>主程序中的读取顺序如下：</p>
<figure class="highlight fortran"><table><tr><td class="gutter"><pre><span class="line">1</span><br><span class="line">2</span><br><span class="line">3</span><br><span class="line">4</span><br><span class="line">5</span><br><span class="line">6</span><br><span class="line">7</span><br><span class="line">8</span><br><span class="line">9</span><br><span class="line">10</span><br><span class="line">11</span><br><span class="line">12</span><br><span class="line">13</span><br><span class="line">14</span><br><span class="line">15</span><br><span class="line">16</span><br><span class="line">17</span><br><span class="line">18</span><br><span class="line">19</span><br><span class="line">20</span><br><span class="line">21</span><br><span class="line">22</span><br><span class="line">23</span><br><span class="line">24</span><br><span class="line">25</span><br><span class="line">26</span><br><span class="line">27</span><br><span class="line">28</span><br><span class="line">29</span><br><span class="line">30</span><br></pre></td><td class="code"><pre><span class="line"> read(iui,*, <span class="keyword">iostat</span>=ierr) isol,inda,indg,imthd,ipol,nda,nds</span><br><span class="line">     <span class="keyword">if</span>(ierr/=<span class="number">0</span>) <span class="keyword">exit</span></span><br><span class="line">     read(iui,*) na0, th0(<span class="number">1</span>:na0)</span><br><span class="line">     read(iui,*) na1,th1(<span class="number">1</span>:na1)</span><br><span class="line">     read(iui,*) nfi,fi(<span class="number">1</span>:nfi)</span><br><span class="line">     read(iui,*) nw0</span><br><span class="line">     read(iui,*) rx(<span class="number">1</span>:nw0+<span class="number">1</span>)</span><br><span class="line">     read(iui,*) rf(<span class="number">1</span>:nw0)</span><br><span class="line">     read(iui,*) nw</span><br><span class="line">     read(iui,*) wl(<span class="number">1</span>:nw)</span><br><span class="line">     read(iui,*) dw(<span class="number">1</span>:nw)</span><br><span class="line">     read(iui,*) galb(<span class="number">1</span>:nw)</span><br><span class="line">     read(iui,*) matm,nlna,nlno</span><br><span class="line">     nlna1=nlna+<span class="number">1</span></span><br><span class="line"><span class="comment">! if ipbf&gt;0 then field out</span></span><br><span class="line">     read(iui,*) ipbf(<span class="number">1</span>:nlna1)</span><br><span class="line"><span class="comment">! input CO2 and ground temperature</span></span><br><span class="line">     read(iui,*) gtmp,co2ppm</span><br><span class="line">     read(iui,*) ifrh,trh</span><br><span class="line">     read(iui,*) npoly,icn,wlcn</span><br><span class="line">     <span class="keyword">do</span> i=<span class="number">1</span>,npoly</span><br><span class="line">        read(iui,*) ncomp(i),cnpt(i)</span><br><span class="line">        read(iui,*) mptc(i,<span class="number">1</span>:ncomp(i))</span><br><span class="line">        read(iui,*) vptc(i,<span class="number">1</span>:ncomp(i))</span><br><span class="line">     <span class="keyword">enddo</span></span><br><span class="line">     <span class="keyword">if</span>(nlno &gt; <span class="number">0</span>) <span class="keyword">then</span></span><br><span class="line">        read(iui,*) dptw(<span class="number">1</span>:nlno)</span><br><span class="line">        read(iui,*) chla(<span class="number">1</span>:nlno)</span><br><span class="line">     <span class="keyword">endif</span></span><br><span class="line"></span><br></pre></td></tr></table></figure>

<p>重复进行其他实验:</p>
<p>说明：</p>
<blockquote>
<ul>
<li>ISOL：0为夜间<br>       1为白天</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>INDA：0为仅通量<br>        1为用于通量和辐射计算</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>INDG：0为朗伯表面<br>        1为用于海洋表面<br>          3为适用于平坦的海洋表面</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>IMTHD：选择通量&#x2F;强度计算模式开关。详见<a target="_blank" rel="noopener external nofollow noreferrer" href="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/%E5%9F%BA%E4%BA%8EFortran90%E7%9A%84Rpstar3%E8%BE%90%E5%B0%84%E4%BC%A0%E8%BE%93%E6%A8%A1%E5%9E%8B%E6%BA%90%E4%BB%A3%E7%A0%81%E6%96%87%E6%A1%A3/MatrixFormulationForTheTransferOfSolarRadiationInAPlane-ParallelScatteringAtmosphere.pdf">Nakajima和Tanaka（1988）Matrix formulation for the transfer of solar radiation in a plane-parallel scattering atmosphere</a></li>
</ul>
</blockquote>
<blockquote>
<p>对于非管理员用户，假定IMTHD&#x3D;2。<br>-1: 计算字段时不进行任何更正。<br>0: DM（Delta-M）-方法。适合通量计算。<br>1: MS（多+单）-方法。适用于辐射计算。<br>2: TMS（截断多个+单个）-方法。有利于反射辐射计算。<br>3: IMS（改进的多+单）-方法。适用于传输辐射计算。<br>（该方法不适用于多孔辐射场）<br>在使用相同的正交流值-NDA-的条件下，与NT&#x2F;DM方法相比，MS&#x2F;TMS&#x2F;IMS方法更精确。<br>在除中厚层以外的大多数情况下，TMS比MS更精确。IMS方法是计算太阳周围强度（太阳光环）的最佳方法。<br>对于极各向异性的相函数，如光学厚度较大的云，TMS将是最好的，且收敛速度比MS方法快得多。</p>
</blockquote>
<p>以下是选择IMTHD的指导表：</p>
<p><img src="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/loading3.gif" data-original="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/%E5%9F%BA%E4%BA%8EFortran90%E7%9A%84Rpstar3%E8%BE%90%E5%B0%84%E4%BC%A0%E8%BE%93%E6%A8%A1%E5%9E%8B%E6%BA%90%E4%BB%A3%E7%A0%81%E6%96%87%E6%A1%A3/imthd.png" alt="IMTHD选择"></p>
<blockquote>
<ul>
<li>IPOL：斯托克斯参数的数量<br>1： 用于标量计算（忽略极化效应）<br>4： 用于矢量计算（计算I、Q、U和V）</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li><p>NDA：半球中高斯正交点的数目，用于积分大气中辐射传输方程的最底角积分。增加该值意味着使用更多计算时间获得更准确的结果（增加NDA**2.5规则）。对计算时间敏感的用户应该通过比较NDA的几个值的结果来调整NDA。</p>
<p>对于通量计算，NDA&#x3D;2-4对于大多数应用来说就足够了。<br>对于辐射率计算，NDA&#x3D;6-10对于TM&#x2F;IMS方法中的大多数应用来说已经足够了。</p>
</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NDS: 半球中用于积分海洋辐射传输方程纳迪亚尔角积分的高斯正交点的数目。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NA0: 太阳天顶角数（大于0）。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>TH0（1:NA0）：太阳天顶角（度）。仅当NA0&gt;0时指定。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NA1：出现辐射的天顶角数（大于0）。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>TH1（1:NA1）：出现辐射的天顶角（度）。<br>透射辐射&lt;90，反射辐射&gt;90。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NFI：方位角数。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>FI（1：NFI）：天顶角（度）。<br>0表示前向散射方向。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NW0：用于计算响应函数的亚波长网格数。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>RX（1：NW0+1）：DW*RX&#x3D;微米波长位移对于子小波网格。<br>-0.5&lt;&#x3D;RX&lt;&#x3D;0.5</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>RF（1:NW0）：辐射计组件的综合响应功能RX功能（相对单位正常）。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>对于单色计算，给出NW0&#x3D;1，RX（1:2）&#x3D;0，RF&#x3D;1。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>平均辐射率：<br>L&#x3D;总和（I&#x3D;1，N；辐射率（I）*射频（I））&#x2F;总和（I&#x3D;1，N；射频（I））</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NWL：中心波长数。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>WL：中心波长（微米）。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>GALB：INDG&#x3D;0时，Lmabert表面的通量反照率。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>如果INDG&#x3D;1，则海面以上10 m处的风速（m&#x2F;s）。<br>（如果u10&gt;0.01，可以进行计算。）</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li><p>MATM：大气数（与AFGL lowtran相同（Kneizys等人，1988年））。</p>
<p>1： 热带，<br>2： 中纬度夏季<br>3：中纬度冬季<br>4： 高纬夏季<br>5：高纬冬季<br>6： 美国标准</p>
</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NLNA：定义大气的子层数。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NLNO：定义海洋的子层数。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>IPBF（1：NLNA+1）：定义要构造的子层的接口号大气。从上到下，即必须按照从上到下的顺序给出：<br>IPBF（1）&#x3D;50，IPBF（NLN+1）&#x3D;1。如果&gt;0，则计算辐射场；如果&lt;0，则无辐射场计算层编号在大气模型中定义MLATMD。</li>
</ul>
</blockquote>
<p>对于IR计算，最好有多个层以获得更高的精度。</p>
<p>默认值如下:</p>
<p>IPBF的表格：</p>
<table>
<thead>
<tr>
<th align="center">接口数量</th>
<th align="center">高度(km)</th>
<th align="center">接口数量</th>
<th align="center">高度(km)</th>
</tr>
</thead>
<tbody><tr>
<td align="center">50</td>
<td align="center">120.0</td>
<td align="center">25</td>
<td align="center">24.0</td>
</tr>
<tr>
<td align="center">49</td>
<td align="center">115.0</td>
<td align="center">24</td>
<td align="center">23.0</td>
</tr>
<tr>
<td align="center">48</td>
<td align="center">110.0</td>
<td align="center">23</td>
<td align="center">22.0</td>
</tr>
<tr>
<td align="center">47</td>
<td align="center">105.0</td>
<td align="center">22</td>
<td align="center">21.0</td>
</tr>
<tr>
<td align="center">46</td>
<td align="center">100.0</td>
<td align="center">21</td>
<td align="center">20.0</td>
</tr>
<tr>
<td align="center">45</td>
<td align="center">95.0</td>
<td align="center">20</td>
<td align="center">19.0</td>
</tr>
<tr>
<td align="center">44</td>
<td align="center">90.0</td>
<td align="center">19</td>
<td align="center">18.0</td>
</tr>
<tr>
<td align="center">43</td>
<td align="center">85.0</td>
<td align="center">18</td>
<td align="center">17.0</td>
</tr>
<tr>
<td align="center">42</td>
<td align="center">80.0</td>
<td align="center">17</td>
<td align="center">16.0</td>
</tr>
<tr>
<td align="center">41</td>
<td align="center">75.0</td>
<td align="center">16</td>
<td align="center">15.0</td>
</tr>
<tr>
<td align="center">40</td>
<td align="center">70.0</td>
<td align="center">15</td>
<td align="center">14.0</td>
</tr>
<tr>
<td align="center">39</td>
<td align="center">65.0</td>
<td align="center">14</td>
<td align="center">13.0</td>
</tr>
<tr>
<td align="center">38</td>
<td align="center">60.0</td>
<td align="center">13</td>
<td align="center">12.0</td>
</tr>
<tr>
<td align="center">37</td>
<td align="center">55.0</td>
<td align="center">12</td>
<td align="center">11.0</td>
</tr>
<tr>
<td align="center">36</td>
<td align="center">50.0</td>
<td align="center">11</td>
<td align="center">10.0</td>
</tr>
<tr>
<td align="center">35</td>
<td align="center">47.5</td>
<td align="center">10</td>
<td align="center">9.0</td>
</tr>
<tr>
<td align="center">34</td>
<td align="center">45.0</td>
<td align="center">9</td>
<td align="center">8.0</td>
</tr>
<tr>
<td align="center">33</td>
<td align="center">42.5</td>
<td align="center">8</td>
<td align="center">7.0</td>
</tr>
<tr>
<td align="center">32</td>
<td align="center">40.0</td>
<td align="center">7</td>
<td align="center">6.0</td>
</tr>
<tr>
<td align="center">31</td>
<td align="center">37.5</td>
<td align="center">6</td>
<td align="center">5.0</td>
</tr>
<tr>
<td align="center">30</td>
<td align="center">35.0</td>
<td align="center">5</td>
<td align="center">4.0</td>
</tr>
<tr>
<td align="center">29</td>
<td align="center">32.5</td>
<td align="center">4</td>
<td align="center">3.0</td>
</tr>
<tr>
<td align="center">28</td>
<td align="center">30.0</td>
<td align="center">3</td>
<td align="center">2.0</td>
</tr>
<tr>
<td align="center">27</td>
<td align="center">27.5</td>
<td align="center">2</td>
<td align="center">1.0</td>
</tr>
<tr>
<td align="center">26</td>
<td align="center">25.0</td>
<td align="center">1</td>
<td align="center">0.0</td>
</tr>
</tbody></table>
<blockquote>
<ul>
<li>IFRH: 0: R.H.从模型大气中使用<br>      1: R.H.使用总R.H.（TRH）</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>TRH: 总相对湿度（0-1）</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NPOLY：粒子多分散体的数量。例如如果有水云和气溶胶，值为2。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li><p>ICN：如果CNPT是每个多分散体的总柱干体积，则为0。</p>
<p>1如果CNPT是每个多分散体的总柱体积<br>2如果CNPT是标度波长WLCN处的光学厚度<br>3如果CNPT是每个波长的光学厚度</p>
</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>WLCN：CNPT作为光学元件的缩放波长（厘米）测量厚度。</li>
</ul>
</blockquote>
<p>重复多分散体的数量</p>
<blockquote>
<ul>
<li>NCOMP：组成多分散体的粒子模型数</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>CNPT：ICN定义的每个多分散体的总体积<br>如果ICN&#x3D;0,1，总体积（cm3&#x2F;cm2）<br>当ICN&#x3D;2,3时，WLCN处的光学厚度</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>MPTC：每个多分散体的粒子型号（请参阅AERDB7一节中的更多信息）</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>VPTC：颗粒模型外部混合的相对干体积浓度*结束重复</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>DPTW：各海洋层的地质厚度【m】</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>叶绿素a：各海洋层的叶绿素a浓度【mg&#x2F;m3】（Bricaud等人（1995），Morel和Maritorena（2001））</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>###粒子模型类型###<br>1： 水<br>2： 冰<br>3： 粉尘类<br>4： 烟尘<br>5： 火山灰<br>6： 75%H2SO4<br>7： 乡村<br>8： 海上浪花<br>9： 城市<br>10：对流层<br>11：黄沙</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>ISPH注释<br>为了使用椭球核数据，需要选择内核数据（FN\u PKRNL1），它取决于球体的纵横比粒子和用户的计算机体系结构（小&#x2F;大端）。内核数据在模块GTPH1（src&#x2F;share&#x2F;01gtph1.f90）中指定，以及可从PKRNL中选择。OUT\u asp{033048069100144207299}in数据库&#x2F;目录。诸如“033”之类的数字表示球状颗粒（“033”表示0.33），由a&#x2F;b定义，其中a是球体长轴的半径，b是短轴的半径。椭球核数据与提供的数据相同rstar7.pack，这是一种RT代码，可从open Slactr网站获得。与rstar7数据的唯一区别是偏振分量。Pstar3除了提供P11数据外，还提供P22、P33、P44、P12和P34数据，配备rstar7。有关内核数据的详细信息，请参见rstar7的描述和Dubovik et al.（2002）的一篇论文。以下read语句用于设置用户定义相对湿度下的水汽分布。<br>读取标题<br>读取（CNG（L，1），L&#x3D;1，NL）</li>
</ul>
</blockquote>
<blockquote>
<p>使用子程序INITS5，在第L层和第M分子的CONG（L，M）数组中设置PPMV单元中的气体浓度分布。Geseus数M定义如下：</p>
<p>1: H2O  2: CO2   3: O3   4: N2O   5: CO   6: CH4   7: O2  8: NO   9: SO2  10: NO2  11: NH3  12: HNO3 13: OH   14: HF  15: HCL  16: HBR  17: HI   18: CLO  19: OCS  20: H2CO  21: HOCL   22: N2  23: HCN  24: CH3CL 25: H2O2  26: C2H2 27: C2H6  28: PH3</p>
<p>在许多应用程序中，用户可能希望在相对湿度，不以PPMV为单位。在主程序中，这种情况假设将相对湿度数据放入数据文件中，然后<br>设定模型大气PPMV浓度后，读取CONG（L，1从数据文件中获取相对湿度值。进入前主程序rstr5，该相对湿度剖面被转换通过常规WVCAL进入PPMV单元。如果你不需要这种情况，只需在主程序中删除此部分。如本例行程序所示，最好修改已装运的干管为您的应用程序编写更好的例程。自结构主程序很简单，用户可以很容易地自己制作一个。</p>
</blockquote>
<h1 id="气体参数文件：MLATMD"><a href="#气体参数文件：MLATMD" class="headerlink" title="气体参数文件：MLATMD"></a>气体参数文件：MLATMD</h1><p>数据文件MLATMD包括常规MLATM读取的Lowtran-7大气模型参数，读取顺序如下：</p>
<figure class="highlight fortran"><table><tr><td class="gutter"><pre><span class="line">1</span><br><span class="line">2</span><br><span class="line">3</span><br><span class="line">4</span><br><span class="line">5</span><br><span class="line">6</span><br><span class="line">7</span><br><span class="line">8</span><br><span class="line">9</span><br><span class="line">10</span><br><span class="line">11</span><br><span class="line">12</span><br><span class="line">13</span><br><span class="line">14</span><br><span class="line">15</span><br><span class="line">16</span><br><span class="line">17</span><br><span class="line">18</span><br><span class="line">19</span><br><span class="line">20</span><br><span class="line">21</span><br><span class="line">22</span><br><span class="line">23</span><br><span class="line">24</span><br><span class="line">25</span><br><span class="line">26</span><br><span class="line">27</span><br><span class="line">28</span><br><span class="line">29</span><br><span class="line">30</span><br></pre></td><td class="code"><pre><span class="line">read NM1,NM2,NATM,NL</span><br><span class="line">    NM=NM1+NM2</span><br><span class="line">    read header</span><br><span class="line">    read AIRM</span><br><span class="line">    read header</span><br><span class="line">    <span class="keyword">DO</span> <span class="number">2</span> I=<span class="number">1</span>,NM</span><br><span class="line">      read IDM(I)</span><br><span class="line">      read NS</span><br><span class="line">      read (IDMS(I,J),J=<span class="number">1</span>,NS)</span><br><span class="line">      read (WMOL(I,J),J=<span class="number">1</span>,NS)</span><br><span class="line">  <span class="number">2</span>   read (RAMS(I,J),J=<span class="number">1</span>,NS)</span><br><span class="line">    read header</span><br><span class="line">    read (ALT(I),I=<span class="number">1</span>,NL)</span><br><span class="line">    <span class="keyword">DO</span> <span class="number">4</span> J=<span class="number">1</span>,NATM</span><br><span class="line">      read header</span><br><span class="line">  <span class="number">4</span>   read (PMATM(I,J),I=<span class="number">1</span>,NL)</span><br><span class="line">    <span class="keyword">DO</span> <span class="number">5</span> J=<span class="number">1</span>,NATM</span><br><span class="line">      read header</span><br><span class="line">  <span class="number">5</span>   read (TMATM(I,J),I=<span class="number">1</span>,NL)</span><br><span class="line">    <span class="keyword">DO</span> <span class="number">6</span> K=<span class="number">1</span>,NM1</span><br><span class="line">    <span class="keyword">DO</span> <span class="number">6</span> J=<span class="number">1</span>,NATM</span><br><span class="line">      read header</span><br><span class="line">  <span class="number">6</span>   read (AMOL(I,K,J),I=<span class="number">1</span>,NL)</span><br><span class="line">    <span class="keyword">DO</span> <span class="number">7</span> J=<span class="number">1</span>,NATM</span><br><span class="line">      read header</span><br><span class="line">  <span class="number">7</span>   read (DNSTY(I,J),I=<span class="number">1</span>,NL)</span><br><span class="line">    <span class="keyword">DO</span> <span class="number">8</span> K=<span class="number">1</span>,NM2</span><br><span class="line">      read header</span><br><span class="line">  <span class="number">8</span>   read (TRAC(I,K),I=<span class="number">1</span>,NL)</span><br><span class="line"></span><br></pre></td></tr></table></figure>



<p>说明：有关详细讨论，请参阅AFGL&#x2F;Lowtran文档</p>
<blockquote>
<ul>
<li>NM1：第一类分子的数目。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NM2：第二类分子的数量。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NATM：模型大气的数量。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NL：高度级别数。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>AIRM：空气分子量</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>IDM（1:NM）：分子代码（A8）。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NS：同位素数量。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>IDMS（1:nm，1:ns）：同位素代码。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>WMOL（1:nm，1:ns）：同位素分子量。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>RAMS（1:nm，1:ns）：同位素相对丰度。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>高度（1:50）：从底部到顶部的高度（km）。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>PMATM（1:nl，1:natm）：水平大气压力（hPa）。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>TMATM（1:nl，1:natm）：水平的大气温度（K）。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>AMOL（1:nl，1:nm1，1:natm）：第一类分子的体积混合比。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>DNSTY（1:nl，1:natm）：模型的空气分子密度（1&#x2F;cm3）。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>TRAC（1:NL，1:NM1）：第二类分子的体积混合比。</li>
</ul>
</blockquote>
<h1 id="粒子参数文件：AERDB"><a href="#粒子参数文件：AERDB" class="headerlink" title="粒子参数文件：AERDB"></a>粒子参数文件：AERDB</h1><p>数据文件AERDB包括由例程GETPAR读取的粒子模型参数，读取顺序如下：</p>
<figure class="highlight fortran"><table><tr><td class="gutter"><pre><span class="line">1</span><br><span class="line">2</span><br><span class="line">3</span><br><span class="line">4</span><br><span class="line">5</span><br><span class="line">6</span><br><span class="line">7</span><br><span class="line">8</span><br><span class="line">9</span><br><span class="line">10</span><br><span class="line">11</span><br><span class="line">12</span><br><span class="line">13</span><br><span class="line">14</span><br><span class="line">15</span><br><span class="line">16</span><br><span class="line">17</span><br><span class="line">18</span><br><span class="line">19</span><br><span class="line">20</span><br><span class="line">21</span><br><span class="line">22</span><br><span class="line">23</span><br><span class="line">24</span><br><span class="line">25</span><br><span class="line">26</span><br><span class="line">27</span><br><span class="line">28</span><br><span class="line">29</span><br><span class="line">30</span><br><span class="line">31</span><br><span class="line">32</span><br><span class="line">33</span><br></pre></td><td class="code"><pre><span class="line">read header</span><br><span class="line">     read RMIN,RMAX</span><br><span class="line">     read NPTC</span><br><span class="line">     <span class="keyword">DO</span> <span class="number">6</span> M=<span class="number">1</span>,NPTC</span><br><span class="line">       read header</span><br><span class="line">   <span class="number">6</span>   read (CNPRF(L,M),L=<span class="number">1</span>,NL)</span><br><span class="line">     <span class="keyword">DO</span> <span class="number">1</span> M=<span class="number">1</span>,NPTC</span><br><span class="line">       read header</span><br><span class="line">       read (ISPCV(I,M),I=<span class="number">1</span>,<span class="number">3</span>)</span><br><span class="line">       read (RFRAC(I,M),I=<span class="number">1</span>,<span class="number">3</span>)</span><br><span class="line">       read RHO(M)</span><br><span class="line">       read NMODE</span><br><span class="line">       DRYAER(<span class="number">1</span>,<span class="number">2</span>,M)=NMODE</span><br><span class="line">       DRYAER(<span class="number">1</span>,<span class="number">3</span>,M)=RMIN</span><br><span class="line">       DRYAER(<span class="number">1</span>,<span class="number">4</span>,M)=RMAX</span><br><span class="line">       <span class="keyword">DO</span> <span class="number">3</span> J=<span class="number">1</span>,NMODE</span><br><span class="line">   <span class="number">3</span>     read (DRYAER(I,J,M),I=<span class="number">2</span>,<span class="number">6</span>)</span><br><span class="line">       read NAW(M)</span><br><span class="line">       <span class="keyword">IF</span>(NAW(M).GT<span class="number">.0</span>) <span class="keyword">then</span></span><br><span class="line">         read header</span><br><span class="line">         <span class="keyword">DO</span> <span class="number">7</span> I=<span class="number">1</span> ,NAW(M)</span><br><span class="line">   <span class="number">7</span>       read (AWCR(I,M,K),K=<span class="number">1</span>,<span class="number">3</span>)</span><br><span class="line">       <span class="keyword">endif</span></span><br><span class="line">   <span class="number">1</span> <span class="keyword">continue</span></span><br><span class="line">     read header</span><br><span class="line">     read NV,NWLV</span><br><span class="line">     read (WLV(I),I=<span class="number">1</span>,NWLV)</span><br><span class="line">     <span class="keyword">DO</span> <span class="number">4</span> IV=<span class="number">1</span>,NV</span><br><span class="line">       read header</span><br><span class="line">       read (RFI(I,IV,<span class="number">1</span>),I=<span class="number">1</span>,NWLV)</span><br><span class="line">   <span class="number">4</span>   read (RFI(I,IV,<span class="number">2</span>),I=<span class="number">1</span>,NWLV)</span><br><span class="line"></span><br><span class="line"></span><br></pre></td></tr></table></figure>



<p>说明：</p>
<blockquote>
<ul>
<li>RMIN：最小粒子半径（cm）。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>RMAX：最大粒子半径（cm）。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>NPTC：粒子模型的数量。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>AERDB中的默认型号集如下所示：</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>粒子模型类型<br>1： 水2：冰<br>3： 粉尘类4：烟尘<br>5： 火山灰6：75%H2SO4<br>7： 乡村8：海上浪花<br>9： 城市10：对流层<br>11： 黄沙</li>
</ul>
</blockquote>
<p>（（CNPRF（L，M），L&#x3D;1，NL），M&#x3D;1，NPTC）：以相对单位表示的干颗粒体积浓度的垂直剖面。<br>（摘自d’Almeida等人）<br>（（ISPCV（I，M），I&#x3D;1,3），M&#x3D;1，NPTC）：我们假设每个模型都有三种成分的内部混合物。ISPCV定义了混合物的基本材料（1-9）。</p>
<blockquote>
<ul>
<li>基础材料（无生长机制）<br>1： 水2：冰<br>3： 粉尘4：水溶性<br>5： 海盐6：烟尘<br>7： 火山灰8：75%H2SO4<br>9： 黄沙</li>
</ul>
</blockquote>
<blockquote>
<p>（（RFRAC（I，M），I&#x3D;1,3），M&#x3D;1，NPTC）：RFRAC定义混合物的干组分体积分数。<br>（ρ（M），M&#x3D;1，NPTC）：干混合物的颗粒密度（g&#x2F;cm3）。<br>（（DRYAER（I，J，M），I&#x3D;2,6，J&#x3D;1，NMODE），M&#x3D;1，NPTC）：卷大小分布（dV&#x2F;dln r）的参数包。<br>其中NMODE是大小分布的模式数。<br>此参数包用于例程VLSPC2定义卷参数包PR的大小分布。</p>
</blockquote>
<blockquote>
<p>PR（1,1）&#x3D;r粒子半径（cm）<br>PR（1,2）&#x3D;模半径的NMODE数<br>PR（1,3）&#x3D;rmin最小粒子半径（cm）<br>PR（1,4）&#x3D;rmax最大粒子半径（cm）<br>对于每个第j个模式（&lt;&#x3D;4）<br>PR（2，j）：模式的功能类型（ITP）。<br>ITP&#x3D;1：幂律<br>PR（3，j）&#x3D;C，PR（4，j）&#x3D;R0，PR（5，j）&#x3D;P<br>vj&#x3D;C<em>（R&#x2F;R0）**（4-P），如果R&gt;R0；&#x3D;C</em>（R&#x2F;R0）</p>
</blockquote>
<blockquote>
<p><strong>4如果R&lt;R0<br>ITP&#x3D;2：对数正常<br>PR（3，j）&#x3D;C，PR（4，j）&#x3D;S，PR（5，j）&#x3D;RM<br>vj&#x3D;C*exp（（ln（R&#x2F;RM）&#x2F;ln（S））</strong>2&#x2F;2）<br>ITP&#x3D;3：修正伽马射线<br>PR（3，j）&#x3D;C，PR（4，j）&#x3D;α，PR（5，j）&#x3D;β，PR（6，j）&#x3D;γ<br>vj&#x3D;C<em>（R1）**（阿尔法+4）经验（-β</em>R1*<em>γ），其中R1&#x3D;R</em>1.0E4<br>（NAW（M），M&#x3D;1，NPTC）：定义颗粒生长的相对湿度数。如果为0，则无增长。<br>（（AWCR（I，M，K），K&#x3D;1,3，I&#x3D;1，NAW（M）），M&#x3D;1，NPTC）：Shettle&amp;Fenn的生长参数（AW、CAW、RMMD）</p>
</blockquote>
<blockquote>
<p>AW水活性解释RMMD的CAW系数。</p>
</blockquote>
<blockquote>
<p>RMMD Hanel的吸水率数据。</p>
</blockquote>
<blockquote>
<p>NV：基本材料的数量。</p>
</blockquote>
<blockquote>
<p>NWLV：定义折射率的波长。</p>
</blockquote>
<blockquote>
<p>（WLV（I），I&#x3D;1，NWLV）：波长，单位为微米。<br>（（RFI（I，IV，1），I&#x3D;1，NWLV），IV&#x3D;1，NV）<br>（（RFI（I，IV，2），I&#x3D;1，NWLV），IV&#x3D;1，NV）：粒子的真实和吸收指数（Mr，Mi）复折射率&#x3D;Mr-i Mi。</p>
</blockquote>
<h1 id="分层设计"><a href="#分层设计" class="headerlink" title="分层设计"></a>分层设计</h1><p>rstar6b采用了三种垂直分层系统。</p>
<blockquote>
<ul>
<li>（1） 分层以定义大气模型。默认为50层layerign公司。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>（2） 分层定义用于辐射传输计算的层系统。可以通过在（1）中的50个层中设置IPBF来设置此分层。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>（3） 如果IPBF&gt;0，则rstar6b将计算通量和强度（如果设置INDA）。</li>
</ul>
</blockquote>
<blockquote>
<p>采用这种相当复杂的分层策略来实现rstar6b在许多应用中的最高效率。一些用户想要计算许多角度方向，而其他一些用户希望<br>设置多个图层。在这些情况下，计算机内存可能会不足如果rstar6b必须记住所有计算值。例如，数组强度的AI0是一个四维数组。对于这些情况，高级用户可以通过调整这些分层来调整最大内存效率保存计算机内存。有关调整方式，请参阅’参数调整’ 。</p>
</blockquote>
<h1 id="主程序结构"><a href="#主程序结构" class="headerlink" title="主程序结构"></a>主程序结构</h1><p>主例程RSTR5需要许多输入数据才能运行。这些数据是读取三个数据文件：data、MLATMD、AERDB和KRNL。退出。设置数据的步骤<br>我们准备了一个数据初始化例程INITS5。</p>
<blockquote>
<ul>
<li>主程序结构如下：<br>（1） 从数据中读取数据<br>（2） 初始化5b<br>（3） pstr4<br>（4） 将结果输出到out</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>用户有两种方法来输入参数，而不是默认值<br>参数：<br>（a） 更改数据文件中的参数。<br>（b） 更换步骤（2）和（3）之间的参数。<br>方法（b）将是暂时更改参数的更简单方法，而方法（a）将是将模型设置为默认值的更好方法参数。</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>主程序RSTR6的输入参数如下：<br>初始化1：初始化例程<br>0：跳过波长独立部分<br>内核文件的IUK设备号<br>INDG-1：无地面<br>0：Lambert曲面<br>1： 初始化的海洋表面<br>2： 无初始化的海面<br>当INDG&gt;0且IMTHD&gt;0时，则进行单次散射校正用于海洋表面反射<br>3： 带初始化的平坦海面<br>4： 无初始化的平坦海面<br>INDA 0：通量<br>1： 通量和辐射<br>IMTHD-1:NT，0:DMS-强度&#x2F;通量法<br>1： MS，2:TMS，3:IMS-强度法。<br>当INDG&gt;0且IMTHD&gt;0时，则进行单次散射校正用于海洋表面反射<br>IMS方法不能应用于IPOL&#x3D;4。<br>ISOL 0：夜间<br>1： 白天<br>NDA：半球中的流数<br>NA0：太阳入射次数<br>TH0（NA0）: 太阳天顶角（0-90度）<br>NA1: 球体中出现的最底角数<br>TH1（NA1）: 紧急最底角（0-90：向下；90-180：向上）NFI方位角数<br>FI（NFI）: 方位角（度）（0方向&#x3D;前向散射）<br>WL: 中心波长（微米）<br>DW: 缩放波长RX的DW缩放因子，单色计算取0<br>NW0： 用于计算响应函数的亚波长网格数<br>RX（NW0）DW*RX&#x3D;亚波长网格的波长位移（微米）-0.5&lt;&#x3D;RX&lt;&#x3D;0.5<br>radiomter的RF（NW0）集成响应功能作为RX的功能（相对单位正常）。<br>对于单色计算，给出NW0&#x3D;1、RX&#x3D;0、RF&#x3D;1。<br>平均辐射率：<br>L&#x3D;总和（I&#x3D;1，N；辐射率（I）*射频（I））&#x2F;总和（I&#x3D;1，N；射频（I））<br>如果INDG&#x3D;0，GALB地面反照率<br>如果INDG&gt;0，则为U10（米&#x2F;秒）<br>ICN 0，1：CNPT为体积<br>2，3:CNPT是WLCN处的光学厚度<br>ICN&#x3D;2，3的WLCN缩放波长<br>NPOLY多分散体数<br>CNPT（NPOLY）总体积（ICN&#x3D;0,1）或光学厚度（ICN&#x3D;2,3）<br>对于每个多分散体。<br>NL定义模型大气的大气层数<br>层间界面高度（km），从下到上<br>各层界面处的PRS（NL）压力（mb）<br>各层界面的TMP（NL）温度（K）<br>NMOL气体数<br>CNG（NL，NMOL）气体浓度（ppmv）<br>CNP（NL，NPOLY）干体积浓度曲线（相对单位）<br>ISPCVP（3，NPOLY）3压缩机基本材料。内部混合物（1-8）<br>RFRACP（3，NPOLY）干混合物的干体积分数<br>ROP（NPOLY）颗粒相对于水的密度<br>干混料的DRYAP（6,4，KPOLY）dV&#x2F;dlnr参数<br>见VLSPC2<br>C值（体积谱系数）是相对的<br>AW的NAWP（NPOLY）编号<br>AWCRP R（NAW，NPOLY，3）1：AW水活性（见Shettle和Fenn。）<br>2： 解释RMMD的CAW系数<br>3： RMMD RMMD<br>基本物种的NV数量（1-8）<br>折射率表的NWLV波长数<br>基波折射率的WLV（NWLV）波长<br>材料（微米）<br>RFI（NWLV，NV，2）基本材料的折射率（mr，mi）<br>&#x3D;mr-i mi<br>具有对数规则波长间隔<br>NLN用于辐射传输的大气层数。<br>IPBF（NLN+1）接口编号，用于定义要构造的子层<br>MLATM中的传输气氛。<br>必须按照从上到下的顺序给出：<br>IPBF（1）&#x3D;50，IPBF（NLN+1）&#x3D;1<br>如果&gt;0，则计算辐射场<br>如果&lt;0，则无辐射场计算</li>
</ul>
</blockquote>
<h1 id="来自R-Pstar3的输出"><a href="#来自R-Pstar3的输出" class="headerlink" title="来自R-Pstar3的输出"></a>来自R-Pstar3的输出</h1><blockquote>
<ul>
<li>THK0（NLN，10）</li>
</ul>
<p>（L，1）：总光学厚度（从大气顶层到地面）<br>（L，2）：粒子光学厚度<br>（L，3）：瑞利光学厚度<br>（L，4）：单次散射反照率</p>
</blockquote>
<blockquote>
<ul>
<li><p>太阳入射辐照度（W&#x2F;m2&#x2F;微米）<br>FLXD0（KNA0，KNTAU）： 向下通量（W&#x2F;m2&#x2F;微米）<br>IPBF&gt;0时</p>
<p>FLXD00（KNA0，KNTAU）:直接入射通量（W&#x2F;m2&#x2F;微米）<br>IPBF&gt;0时</p>
<p>FLXU0（KNA0，KNTAU）:向上通量（W&#x2F;m2&#x2F;微米）<br>IPBF&gt;0时</p>
<p>AI0（KNA1U、KNA0、KNFI、KNTAU、1）<br>辐射率:（W&#x2F;m2&#x2F;微米&#x2F;str）<br>IPBF&gt;0时</p>
</li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>PL:（KNA1U、KNA0、KNFI、KNTAU）<br>偏振度&#x3D;sqrt（Q<strong>2+U</strong>2+V**2）&#x2F;I</li>
</ul>
</blockquote>
<h1 id="调整参数大小"><a href="#调整参数大小" class="headerlink" title="调整参数大小"></a>调整参数大小</h1><blockquote>
<ul>
<li>您可能需要更改应用程序的参数设置。<br>参数的含义如下：<br>KNA0：最大太阳入射次数。<br>KNA1U：球体中的最大出射角数。<br>KNFI：最大方位角数。<br>KNW0：最大接收数。<br>KNW：最大波长数。<br>KNLN：辐射传输计算的最大层数。它应该等于或小于49。<br>KNTAU：辐射场所在层界面的最大数量计算它应等于或小于KNLN。<br>以下参数调整需要详细了解rstar6b编码,不建议普通用户更改。<br>KNDM：辐射传输的最大流数。<br>KNL：模型大气水平的最大数量（一般用户不应更改此设置）<br>KNM0：最大GESE数。<br>KPOLY：辐射传输的最大粒子多色散数。<br>KAW：水活动的最大数量。<br>KWLV：定义折射率的最大波长数<br>基本材料。<br>KNV：粒子基本材料的最大数量。<br>KLGN1：相位函数勒让德展开的最大阶数+1<br>在IMS方法中。<br>kptc：粒子类型的最大数量<br>kp：用于吸收的压力格栅数量（与ckd.g\U ch2相关）<br>kt：用于吸收的温度网格数（与ckd.g\U ch2相关）<br>kmol：吸收气体种类数（与ckd、g\U ch2相关）<br>knang：散射角网格数（与KRNL.OUT、pyphsf、dkrnl相关）<br>kplk1：多项式拟合普朗克函数的最大阶数<br>kww：气体吸收的最大波数积分数<br>kch：k分布的积分点数（与ckd.g\U ch2相关）<br>kintvl：KRNL的最大大小参数网格数。外出<br>kpol：最大极化数（&#x3D;1表示rstar）<br>krvis:pyphsf\u vis的波长网格数<br>krir:pyphsf\u ir的波长网格数<br>kintvl\U du：dkrnl的尺寸参数网格数<br>knm：最大气体种类数</li>
</ul>
</blockquote>
<h1 id="注意事项"><a href="#注意事项" class="headerlink" title="注意事项"></a>注意事项</h1><blockquote>
<ul>
<li>此包是从Rstar7 ocean修改而来的，因此存在一些差异来自下面列出的原始Pstar3包。<br>*LBL计算部分不包含此包。<br>*IMS方法不适用于极化辐射场。<br>*太阳光谱仅包含2个（将更新）。<br>*平坦的海面对于蒸发辐射场来说是不好的。<br>（可能有一些未知的bug…）</li>
</ul>
</blockquote>
<h1 id="参考文献"><a href="#参考文献" class="headerlink" title="参考文献"></a>参考文献</h1><p>[1]d’Almeida, G. A., P. Koepke, and E. P. Shettle, 1991: Atmospheric aerosols.</p>
<p> Global climatology and radiative characteristics. A. Deepak Publishing.</p>
<p>[2]Dubovik, O., B. N. Holben, T. Lapyonok, A. Sinyuk, M. I. Mishchenko, P. Yang,</p>
<p> and I. Slutsker 2002: Non-spherical aerosol retrieval method employing </p>
<p> light scattering by spheroids, Geophys. Res. Lett., 29(10), 1415, </p>
<p> doi:10.1029&#x2F;2001GL014506.</p>
<p>[3]Hanel, G., 1976: The properties of atmospheric aerosol particles as functions</p>
<p> of the relative humidity as a thermodynamic equillibrium with the surrounding</p>
<p> moist air. Advances in Geophys., 19, 73-188.</p>
<p>[4]Kneizys, F. X., E. P. Shettle, L. W. Abreu, J. H. Chetwynd, G. P. Anderson,</p>
<p> W. O. Gallery, J. E. A. Selby, and S. A. Clough, 1988: Users Guide to</p>
<p> LOWTRAN 7. AFGL-TR-88-0177.</p>
<p>[5]Nakajima, T., and M. Tanaka, 1986: Matrix formulations for the transfer</p>
<p> of solar radiation in a plane-parallel scattering atmosphere. J.</p>
<p> Quant. Spectrosc. Radiat. Transfer, 35, 13-21.</p>
<p>[6]Nakajima, T., and M. Tanaka, 1988: Algorithms for radiative intensity</p>
<p> calculations in moderately thick atmospheres using a truncation</p>
<p> approximation. J. Quant. Spectrosc. Radiat. Transfer , 40, 51-69.</p>
<p>[7]Nakajima, T., M. Tsukamoto, Y. Tsushima, A. Numaguti, and T. Kimura, 2000: </p>
<p> Modeling of the radiative process in an atmospheric general circulation model.</p>
<p> Appl. Opt., 39, 4869-4878.</p>
<p>[8]Ota, Y., A. Higurashi, T. Nakajima and T. Yokota, 2010: Matrix formulations</p>
<p> of radiative transfer including the polarization effect in a coupled</p>
<p> atmosphere-ocean system, J. Quant. Spectrosc. Radiat. Transfer,</p>
<p> 111, 878-894, doi:10.1016&#x2F;j.jqsrt.2009.11.021.</p>
<p>[9]Sekiguchi, M., and T. Nakaima, 2008: A study of the absorption process and</p>
<p> its computational optimization in an atmospheric general circulation model.</p>
<p> J. Quant. Spectrosc. Radiat. Transfer, 109, 2779-2793,</p>
<p> doi:10.1016&#x2F;j.jqsrt.2008.07.013.</p>
<p>[10]Shettle, E. P., and R. W. Fenn, 1979: Models for the aerosols of the lower</p>
<p> atmosphere and the effects of humidity variations on their optical</p>
<p> properties. AFGL-TR-79-0214.</p>
<p>[11]Yang, P., K. N. Liou, K. Wyser, and D. Mitchell 2000: Parameterization of </p>
<p> the scattering and absorption properties of individual ice crystals, </p>
<p> J. Geophys. Res., 105(D4), 4699–4718, doi:10.1029&#x2F;1999JD900755.</p>
<p>[12]Yang, P., H. Wei, H.-L. Huang, B. A. Baum, Y. X. Hu, G. W. Kattawar, </p>
<p> M. I. Mishchenko, and Q. Fu, 2005: Scattering and absorption property </p>
<p> database for nonspherical ice particles in the near- through far-infrared</p>
<p> spectral region Appl. Opt. 44, 5512-5523.</p>
<h1 id="出资者"><a href="#出资者" class="headerlink" title="出资者"></a>出资者</h1><blockquote>
<ul>
<li>Higurashi，Akiko（海面）<br>国家环境研究所<br>电话：+81-29-850-2423；传真：+81-29-850-2960<br>电子邮件：<a href="mailto:&#104;&#x61;&#x6b;&#x69;&#107;&#x6f;&#64;&#110;&#105;&#x65;&#115;&#x2e;&#103;&#x6f;&#x2e;&#x6a;&#112;" rel="external nofollow noreferrer">&#104;&#x61;&#x6b;&#x69;&#107;&#x6f;&#64;&#110;&#105;&#x65;&#115;&#x2e;&#103;&#x6f;&#x2e;&#x6a;&#112;</a></li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>中岛，Teruyuki（辐射传输）<br>国家环境研究所<br>电子邮件：<a href="mailto:&#116;&#x65;&#114;&#114;&#121;&#45;&#x6e;&#107;&#x6a;&#64;&#110;&#x69;&#102;&#116;&#121;&#46;&#x63;&#x6f;&#x6d;" rel="external nofollow noreferrer">&#116;&#x65;&#114;&#114;&#121;&#45;&#x6e;&#107;&#x6a;&#64;&#110;&#x69;&#102;&#116;&#121;&#46;&#x63;&#x6f;&#x6d;</a></li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>Ota，Yoshifumi（海洋辐射传输）<br>日本气象厅<br>电子邮件：<a href="mailto:&#111;&#116;&#97;&#121;&#111;&#115;&#x68;&#105;&#64;&#109;&#x65;&#x74;&#x2e;&#x6b;&#105;&#x73;&#104;&#111;&#117;&#46;&#103;&#111;&#x2e;&#x6a;&#112;" rel="external nofollow noreferrer">&#111;&#116;&#97;&#121;&#111;&#115;&#x68;&#105;&#64;&#109;&#x65;&#x74;&#x2e;&#x6b;&#105;&#x73;&#104;&#111;&#117;&#46;&#103;&#111;&#x2e;&#x6a;&#112;</a></li>
</ul>
</blockquote>
<blockquote>
<ul>
<li>Sekiguchi，Miho（气体吸收模型和辐射传输）<br>东京海洋科技大学<br>电话和传真：+81-3-5245-7463<br>电子邮件：<a href="mailto:&#109;&#x69;&#x68;&#x6f;&#64;&#x6b;&#x61;&#105;&#x79;&#x6f;&#x64;&#97;&#x69;&#x2e;&#97;&#99;&#46;&#106;&#112;" rel="external nofollow noreferrer">&#109;&#x69;&#x68;&#x6f;&#64;&#x6b;&#x61;&#105;&#x79;&#x6f;&#x64;&#97;&#x69;&#x2e;&#97;&#99;&#46;&#106;&#112;</a></li>
</ul>
</blockquote>
</article><div class="post-copyright"><div class="post-copyright__author"><span class="post-copyright-meta">文章作者: </span><span class="post-copyright-info"><a href="mailto:guojiaxiang0820@gmail.com" rel="external nofollow noreferrer">洛沐</a></span></div><div class="post-copyright__type"><span class="post-copyright-meta">文章链接: </span><span class="post-copyright-info"><a href="https://www.guojxblog.cn/archives/8c1f6a98.html">https://www.guojxblog.cn/archives/8c1f6a98.html</a></span></div><div class="post-copyright__notice"><span class="post-copyright-meta">版权声明: </span><span class="post-copyright-info">本博客所有文章除特别声明外，均采用 <a href="https://creativecommons.org/licenses/by-nc-sa/4.0/" rel="external nofollow noreferrer" target="_blank">CC BY-NC-SA 4.0</a> 许可协议。转载请注明来自 <a href="https://www.guojxblog.cn" target="_blank">洛沐の人间客栈</a>！</span></div></div><div class="tag_share"><div class="post-meta__tag-list"><a class="post-meta__tags" href="/tags/%E8%BE%90%E5%B0%84%E4%BC%A0%E8%BE%93/">辐射传输</a></div><div class="post_share"><div class="social-share" data-image="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/%E5%9F%BA%E4%BA%8EFortran90%E7%9A%84Rpstar3%E8%BE%90%E5%B0%84%E4%BC%A0%E8%BE%93%E6%A8%A1%E5%9E%8B%E6%BA%90%E4%BB%A3%E7%A0%81%E6%96%87%E6%A1%A3/rpstar.jpeg" data-sites="facebook,twitter,wechat,weibo,qq"></div><link rel="stylesheet" href="https://cdn.jsdelivr.net/npm/social-share.js/dist/css/share.min.css" media="print" onload="this.media='all'"><script src="https://cdn.jsdelivr.net/npm/social-share.js/dist/js/social-share.min.js" defer></script></div></div><div class="post-reward"><div class="reward-button"><i class="fas fa-qrcode"></i> 打赏</div><div class="reward-main"><ul class="reward-all"><li class="reward-item"><a href="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/wechat.jpeg" rel="external nofollow noreferrer" target="_blank"><img class="post-qr-code-img" src="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/loading3.gif" data-original="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/wechat.jpeg" alt="微信"/></a><div class="post-qr-code-desc">微信</div></li><li class="reward-item"><a href="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/alipay.jpeg" rel="external nofollow noreferrer" target="_blank"><img class="post-qr-code-img" src="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/loading3.gif" data-original="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/alipay.jpeg" alt="支付宝"/></a><div class="post-qr-code-desc">支付宝</div></li></ul></div></div><nav class="pagination-post" id="pagination"><div class="prev-post pull-left"><a href="/archives/c144f20a.html"><img class="prev-cover" src="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/loading3.gif" data-original="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/bg.jpeg" onerror="onerror=null;src='/img/404.jpg'" alt="cover of previous post"><div class="pagination-info"><div class="label">上一篇</div><div class="prev_info">敢于重新开始</div></div></a></div><div class="next-post pull-right"><a href="/archives/49977e84.html"><img class="next-cover" src="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/loading3.gif" data-original="https://luomublog.oss-cn-qingdao.aliyuncs.com/ImgHost/beauty_top.jpeg" onerror="onerror=null;src='/img/404.jpg'" alt="cover of next post"><div class="pagination-info"><div class="label">下一篇</div><div class="next_info">惊人的美</div></div></a></div></nav><hr/><div id="post-comment"><div class="comment-head"><div class="comment-headline"><i class="fas fa-comments fa-fw"></i><span> 评论</span></div></div><div class="comment-wrap"><div><div id="lv-container" data-id="city" data-uid="MTAyMC81NjIzOS8zMjcwMg=="></div></div></div></div></div><div class="aside-content" id="aside-content"><div class="sticky_layout"><div class="card-widget" id="card-toc"><div class="item-headline"><i class="fas fa-stream"></i><span>目录</span><span class="toc-percentage"></span></div><div class="toc-content is-expand"><ol class="toc"><li class="toc-item toc-level-1"><a class="toc-link" href="#%E2%80%9CR-Pstar3%E2%80%9D%E7%AE%80%E4%BB%8B%EF%BC%9F"><span class="toc-number">1.</span> <span class="toc-text">“R-Pstar3”简介？</span></a></li><li class="toc-item toc-level-1"><a class="toc-link" href="#%E8%BD%AF%E4%BB%B6%E5%8C%85%E7%A8%8B%E5%BA%8F%E6%B8%85%E5%8D%95"><span class="toc-number">2.</span> <span class="toc-text">软件包程序清单</span></a></li><li class="toc-item toc-level-1"><a class="toc-link" href="#%E6%BA%90%E7%A0%81%E8%BF%90%E8%A1%8C%E4%B8%8E%E7%BC%96%E8%AF%91"><span class="toc-number">3.</span> <span class="toc-text">源码运行与编译</span></a></li><li class="toc-item toc-level-1"><a class="toc-link" href="#%E5%AE%9E%E9%AA%8C%E5%8F%82%E6%95%B0%E6%96%87%E4%BB%B6%EF%BC%9A%E6%95%B0%E6%8D%AE"><span class="toc-number">4.</span> <span class="toc-text">实验参数文件：数据</span></a></li><li class="toc-item toc-level-1"><a class="toc-link" href="#%E6%B0%94%E4%BD%93%E5%8F%82%E6%95%B0%E6%96%87%E4%BB%B6%EF%BC%9AMLATMD"><span class="toc-number">5.</span> <span class="toc-text">气体参数文件：MLATMD</span></a></li><li class="toc-item toc-level-1"><a class="toc-link" href="#%E7%B2%92%E5%AD%90%E5%8F%82%E6%95%B0%E6%96%87%E4%BB%B6%EF%BC%9AAERDB"><span class="toc-number">6.</span> <span class="toc-text">粒子参数文件：AERDB</span></a></li><li class="toc-item toc-level-1"><a class="toc-link" href="#%E5%88%86%E5%B1%82%E8%AE%BE%E8%AE%A1"><span class="toc-number">7.</span> <span class="toc-text">分层设计</span></a></li><li class="toc-item toc-level-1"><a class="toc-link" href="#%E4%B8%BB%E7%A8%8B%E5%BA%8F%E7%BB%93%E6%9E%84"><span class="toc-number">8.</span> <span class="toc-text">主程序结构</span></a></li><li class="toc-item toc-level-1"><a class="toc-link" href="#%E6%9D%A5%E8%87%AAR-Pstar3%E7%9A%84%E8%BE%93%E5%87%BA"><span class="toc-number">9.</span> <span class="toc-text">来自R-Pstar3的输出</span></a></li><li class="toc-item toc-level-1"><a class="toc-link" href="#%E8%B0%83%E6%95%B4%E5%8F%82%E6%95%B0%E5%A4%A7%E5%B0%8F"><span class="toc-number">10.</span> <span class="toc-text">调整参数大小</span></a></li><li class="toc-item toc-level-1"><a class="toc-link" href="#%E6%B3%A8%E6%84%8F%E4%BA%8B%E9%A1%B9"><span class="toc-number">11.</span> <span class="toc-text">注意事项</span></a></li><li class="toc-item toc-level-1"><a class="toc-link" href="#%E5%8F%82%E8%80%83%E6%96%87%E7%8C%AE"><span class="toc-number">12.</span> <span class="toc-text">参考文献</span></a></li><li class="toc-item toc-level-1"><a class="toc-link" href="#%E5%87%BA%E8%B5%84%E8%80%85"><span class="toc-number">13.</span> <span class="toc-text">出资者</span></a></li></ol></div></div></div></div></main><footer id="footer"><div id="footer-wrap"><div class="copyright">&copy;2021 - 2023  <i id="heartbeat" class="fa fas fa-heartbeat"></i> 洛沐</div><div class="footer_custom_text">谢谢你来看<a href="https://www.guojxblog.cn/" style='color:red;Font-size:36'>我</a>，你今天真好看😘</div></div><link rel="stylesheet" href="https://cdn.jsdelivr.net/gh/HCLonely/images@master/others/heartbeat.min.css"></footer></div><div id="rightside"><div id="rightside-config-hide"><button id="readmode" type="button" title="阅读模式"><i class="fas fa-book-open"></i></button><button id="translateLink" type="button" title="简繁转换">簡</button><button id="darkmode" type="button" title="浅色和深色模式转换"><i class="fas fa-adjust"></i></button><button id="hide-aside-btn" type="button" title="单栏和双栏切换"><i class="fas fa-arrows-alt-h"></i></button></div><div id="rightside-config-show"><button id="rightside_config" type="button" title="设置"><i class="fas fa-cog fa-spin"></i></button><button class="close" id="mobile-toc-button" type="button" title="目录"><i class="fas fa-list-ul"></i></button><a id="to_comment" href="#post-comment" title="直达评论"><i class="fas fa-comments"></i></a><button id="go-up" type="button" title="回到顶部"><i class="fas fa-arrow-up"></i></button></div></div><div id="local-search"><div class="search-dialog"><nav class="search-nav"><span class="search-dialog-title">本地搜索</span><span id="loading-status"></span><button class="search-close-button"><i class="fas fa-times"></i></button></nav><div class="is-center" id="loading-database"><i class="fas fa-spinner fa-pulse"></i><span>  数据库加载中</span></div><div class="search-wrap"><div id="local-search-input"><div class="local-search-box"><input class="local-search-box--input" placeholder="搜索文章" type="text"/></div></div><hr/><div id="local-search-results"></div></div></div><div id="search-mask"></div></div><div><script src="/js/utils.js"></script><script src="/js/main.js"></script><script src="/js/tw_cn.js"></script><script src="https://cdn.jsdelivr.net/npm/medium-zoom/dist/medium-zoom.min.js"></script><script src="https://cdn.jsdelivr.net/npm/instant.page/instantpage.min.js" type="module"></script><script src="https://cdn.jsdelivr.net/npm/node-snackbar/dist/snackbar.min.js"></script><script>function panguFn () {
  if (typeof pangu === 'object') pangu.autoSpacingPage()
  else {
    getScript('https://cdn.jsdelivr.net/npm/pangu/dist/browser/pangu.min.js')
      .then(() => {
        pangu.autoSpacingPage()
      })
  }
}

function panguInit () {
  if (false){
    GLOBAL_CONFIG_SITE.isPost && panguFn()
  } else {
    panguFn()
  }
}

document.addEventListener('DOMContentLoaded', panguInit)</script><script src="/js/search/local-search.js"></script><div class="js-pjax"><script>function loadLivere () {
  if (typeof LivereTower === 'object') {
    window.LivereTower.init()
  }
  else {
    (function(d, s) {
        var j, e = d.getElementsByTagName(s)[0];
        if (typeof LivereTower === 'function') { return; }
        j = d.createElement(s);
        j.src = 'https://cdn-city.livere.com/js/embed.dist.js';
        j.async = true;
        e.parentNode.insertBefore(j, e);
    })(document, 'script');
  }
}

if ('Livere' === 'Livere' || !true) {
  if (true) btf.loadComment(document.getElementById('lv-container'), loadLivere)
  else loadLivere()
}
else {
  function loadOtherComment () {
    loadLivere()
  }
}</script></div><script src="/js/mobile_side.js"></script><script id="click-heart" src="https://cdn.jsdelivr.net/npm/butterfly-extsrc@1/dist/click-heart.min.js" async="async" mobile="true"></script><div class="pjax-reload"><script async="async">var arr = document.getElementsByClassName('recent-post-item');
for(var i = 0;i<arr.length;i++){
    arr[i].classList.add('wow');
    arr[i].classList.add('animate__zoomIn');
    arr[i].setAttribute('data-wow-duration', '1s');
    arr[i].setAttribute('data-wow-delay', '700ms');
    arr[i].setAttribute('data-wow-offset', '100');
    arr[i].setAttribute('data-wow-iteration', '1');
}</script><script async="async">var arr = document.getElementsByClassName('card-widget');
for(var i = 0;i<arr.length;i++){
    arr[i].classList.add('wow');
    arr[i].classList.add('animate__zoomIn');
    arr[i].setAttribute('data-wow-duration', '');
    arr[i].setAttribute('data-wow-delay', '');
    arr[i].setAttribute('data-wow-offset', '');
    arr[i].setAttribute('data-wow-iteration', '');
}</script></div><script defer="defer" src="https://cdn.jsdelivr.net/gh/graingert/wow@1.3.0/dist/wow.min.js"></script><script defer="defer" src="/js/wow_init.js"></script><script async data-pjax src="//busuanzi.ibruce.info/busuanzi/2.3/busuanzi.pure.mini.js"></script></div><script>
            window.imageLazyLoadSetting = {
                isSPA: true,
                preloadRatio: 3,
                processImages: null,
            };
        </script><script>window.addEventListener("load",function(){var t=/\.(gif|jpg|jpeg|tiff|png)$/i,r=/^data:image\/[a-z]+;base64,/;Array.prototype.slice.call(document.querySelectorAll("img[data-original]")).forEach(function(a){var e=a.parentNode;"A"===e.tagName&&(e.href.match(t)||e.href.match(r))&&(e.href=a.dataset.original)})});</script><script>!function(n){n.imageLazyLoadSetting.processImages=o;var e=n.imageLazyLoadSetting.isSPA,i=n.imageLazyLoadSetting.preloadRatio||1,r=Array.prototype.slice.call(document.querySelectorAll("img[data-original]"));function o(){e&&(r=Array.prototype.slice.call(document.querySelectorAll("img[data-original]")));for(var t,a=0;a<r.length;a++)0<=(t=(t=r[a]).getBoundingClientRect()).bottom&&0<=t.left&&t.top<=(n.innerHeight*i||document.documentElement.clientHeight*i)&&function(){var t,e,n,i,o=r[a];t=o,e=function(){r=r.filter(function(t){return o!==t})},n=new Image,i=t.getAttribute("data-original"),n.onload=function(){t.src=i,e&&e()},t.src!==i&&(n.src=i)}()}o(),n.addEventListener("scroll",function(){var t,e;t=o,e=n,clearTimeout(t.tId),t.tId=setTimeout(function(){t.call(e)},500)})}(this);</script><script async>window.onload=function(){var a=document.createElement('script'),b=document.getElementsByTagName('script')[0];a.type='text/javascript',a.async=!0,a.src='/sw-register.js?v='+Date.now(),b.parentNode.insertBefore(a,b)};</script></body></html>